Deconfinement transition and Luttinger to Fermi Liquid crossover in quasi 

one-dimensional systems 
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We investigate a system of one dimensional Hubbard chains of interacting fermions coupled by 
inter-chain hopping. Using a generalization of the Dynamical Mean Field Theory we study the 
deconfinement transition from a Mott insulator to a metal and the crossover between Luttinger 
and Fermi liquid phases. One-particle properties, local spin response and inter-chain optical 
conductivity are calculated. Possible applications to organic conductors are discussed. 
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The nature of the metallic phase of interacting electron 
systems depends strongly on dimensionality. In three di- 
mensions, Fermi liquid (FL) theory applies, whereas in 
one dimension the quasi-particle concept breaks down, 
leading to a different kind of low-energy fixed point 
known as a Luttinger liquid (LL). For commensurate 
electron fillings, strong enough repulsive interactions de- 
stroy the metallic state altogether by opening a Mott gap. 
This phenomenon exists in all dimensions but the one- 
dimensional case is particularly favorable |l|. In quasi 
one-dimensional (QID) systems, interchain hopping can 
induce a (deconfinement) transition from the Mott insu- 
lating (MI) state to a metallic state and crossovers be- 
tween different metallic behaviors (Fig. 1). Besides their 
intrinsic theoretical interest, understanding these phe- 
nomena is directly relevant for a number of compounds 
such as the organic (super)-conductors (Bechgaard salts), 
which are three dimensional stacks of quarter- filled 
chains [§• Indeed, in these compounds some of the low 
temperature properties are well described by Fermi liq- 
uid theory, whereas optical [^ and transport properties 
[^ have shown that the high temperature phase is either 
a Luttinger liquid or a Mott insulator. 

Describing these QID systems is not an easy task. The 
transverse hopping t± is a relevant perturbation on the 
LL 1^, I3, |7|. Hence, perturbative renormalization group 
calculations yield an estimate of the crossover scale [g], ||] 
but fail below that scale. In the MI state, electrons are 
confined on the chains by the Mott gap. A finite criti- 
cal value of i_L is needed to induce an insulator to metal 
transition. Such a transition has been advocated ^, 10 




to 

be at the root of the change from insulating to metal- 
lic behavior observed in Bechgaard salts when increasing 
pressure or when going from the TMTTF to the TMTSF 
family. Thus the effects of the inter-chain hopping that 
are the most important physically cannot be handled 
reliably by perturbative methods. Although some non- 



FIG. 1: Schematic phase diagram of a quasi one-dimensional 
system, as a function of temperature T and inter-chain hop- 
ping t±. The Mott-Insulator (MI), Luttinger Liquid (LL) 
and Fermi Liquid (FL) regimes are displayed. All lines are 
crossovers, except for the T = transition from a MI to a 
FL. The MI phase is present only for commensurate fillings. 
Possible phases with long-range order have been omitted. 



perturbative studies of deconfinement have been made 
for a finite number of chains [[ll| the case of an infinite 
system is still open. Some of the key questions yet to be 
answered are: (i) What is the crossover scale from the 
LL to MI, or from the LL to FL (Fig.l), and is there 
only one crossover scale for the different physical proper- 
ties (transport, spin response, single-particle properties 
etc..) ? (ii) What is the nature of the low-temperature 
FL state, and in particular is the shape of the Fermi 
surface (FS) affected by interactions ? Do quasi-particle 
(QP) properties such as the quasi-particle residue Z{k), 
vary significantly as the wavevector k moves along the 
FS ? 

In order to deal with these questions in a non- 
perturbative manner, a new method has been pro- 



posed jij, |l^ . This method, which generahzes the single 
site dynamical mean- field theory (DMFT) [|4|, replaces 
the quasi one-dimensional system by a single effective 
chain from which electrons can hop to a self-consistent 
bath. We thus nickname it chain-DMFT (ch-DMFT). 
It is exact in the limit of decoupled chains, while it 
reduces to usual DMFT in the opposite limit of zero- 
hopping along the chains, hence capturing the physics 
of a strongly correlated FL (or MI) WM. Solving the 
mean field equations is a very challenging task however. 
In this paper, the first quantitative solution of the ch- 
DMFT equations is presented. We show that the above 
issues can be valuably addressed, and some of the above 
questions answered. In view of the difficulty of the prob- 
lem only the Hubbard model is considered in this paper, 
but more realistic extensions are in sight, as discussed 
below. 

Let us consider a system of coupled chains described 
by a sum of Hamiltonians of the isolated chains plus 
a hopping term between neighboring chains. The ch- 
DMFT approximates this system by the effective action: 
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where -ff JJ) i^ ^^^ interacting part of the on-chain Hamil- 
tonian. Qq is an effective propagator determined from a 
self-consistency condition: the Green's function G{i — 
j,T — t') = —{c{i,T)c'^{j,T'))cS calculated from ScS 
should coincide with the on-chain Green's function of 
the original problem, with the same self-energy S = 
Q^^ -G-^. This reads: 



and we consider a two-dimensional array of chains form- 
ing a square lattice with nearest neighbor hopping, so 
that the transverse dispersion is e^(fc^) = —2t_\_ cosfcj,. 
We solve the effective chain problem for chains of 
16 or 32 sites with periodic boundary conditions using 
the QMC Hirsch-Fye algorithm ||l|l. Using 32 time- 
slices in imaginary time allows to reliably access tem- 
peratures down to T/W ~ 1/50 = 0.02 with U/W ~ 1 
where W — 41 is the bandwidth of the one-dimensional 
chain. Typically, about 5000 QMC sweeps and 10 to 
15 ch-DMFT iterations are sufficient to reach conver- 
gence. The main quantities that we calculate and ana- 
lyze are: a) On-chain (i.e summed over fcj^) single particle 
Green's functions: G'(fc,zw„), and self-energies I](fc,iw„) 
(obtained from Q^^ — G^^ [21])- This allows to iden- 
tify the location of the FS k^{k) by solving e±{k^) — 
jjL — ReS(fc, za;„=i) — e-k and the QP residue Z = Z{k^) 
by fitting the slope (= 1 — Z^^) of S vs. iuJn at a spe- 
cific FS point, b) On-chain spin and charge response 
functions and in particular the local spin response func- 
tion: xs{r) = {SHj,0)SHj,t)) = Ek,k±Xs{k,k^,T). 
Indeed, in a LL (or with Kp = 1 in a FL) we have 
XsiT) — Xs(/3/2) (siuTrr//?)^^ ^' in the asymptotic 
regime where (3, t are larger than some coherence scale. 
By fitting our results to this form, we can define an ef- 
(1) fective Kp. Furthermore, the temperature dependence 
of Xs{P/'2) is closely related to the NMR relaxation rate 
l/Ti = lim^^o x"('^)/'-^- Ii^ a FL liquid at low-enough 
temperature @: 1/Ti = /3xs(/3/2)/27r2. In a LL 1/Ti 
and /3xs(/3/2) have the same T-dependence (albeit with 
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' iuJn + l-J.- Cfc - S(fc, iuJn) - e±_ 



(2) 



Let €_\_{k±) denote the Fourier transform of the inter- 
chain hopping t^^', D{e±) = E^a ^[^^ - ^^ik±)] 
the corresponding density of states, k the momen- 
tum in the chain direction and LOn the Matsubara 
frequencies corresponding to the inverse temperature 
/3. The ch-DMFT equations (0j|) fully determine the 
self-energy and Green's function of the coupled chains 



G{k,k_L,iujn)~ 



+ pL ~ ek ~ e_L(fc_L) - T.{k,iu>n)- 



The ch-DMFT approach can be rigorously justified in the 
formal limit where the lattice connectivity in the trans- 
verse direction is taken to infinity. In systems with finite 
transverse connectivity it can be viewed as an approxi- 
mation neglecting the dependence of the self-energy on 
transverse momentum, keeping both frequency and in- 
chain momentum dependence: S = S(fc,ia;„). In the 
following, we specialize to the Hubbard model, i.e for 
each chain: 
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chain optical conductivity. Within ch-DMFT vertex cor- 
rections drop out ||l^ and the conductivity is given by 

a±{iu;) ex -^ — — sin^ A;^ / 7;- :^ T^ G{k,k±,iu;„i) 

LU I ZTT I 2tT P ^-^ 

m 

xG{k, k±,iLu„i + iuj) (4) 

Note that we have taken into account the fc_L- dependence 
of the current vertex in (0). In practice, we perform an- 
alytically the fc^-integration in (0) and numerically the 
k-integration and Matsubara summation. We then use 
the Maximum Entropy (MaxEnt) algorithm |I^ to con- 
tinue a^{iuj) to the real axis and obtain the inter-chain 
optical conductivity Re(Tx((jj,T). 

We first discuss our results away from half-filling, for 
U/W — 1 and pL = 0.2 (corresponding to a total density 
n ~ 0.8). The numerical method was thoroughly tested 
for decoupled chains (tj^ = 0), for which our fitting pro- 
cedure of Xs yields Kp ~ 0.7, in agreement with the exact 
result |1^ . We then study how inter-chain coherence de- 
velops as temperature is lowered for coupled chains with 
t^/W = 0.14. Fig. H displays the effective Kp as a func- 
tion of temperature: a clear crossover from a LL (with 
Kp ~ 0.7) to a FL [Kp = 1) is seen as temperature is 
lowered below T* /W ~ 1/44. According to ||, §, the 
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FIG. 2: Effective Kp vs. temperature in the doped case 
(filling n ~ 0.8) for U/W = 1.0, t±/W = 0.14. 



TABLE I: QP weights Z{kj_) for different points on the FS 
(doped case: n ~ 0.8, U/W = 1.0, t^/W = 0.14). 

fcx/TT 0.19 0.36 0.50 0.62 0.74 0.91 
Z{k±_) 0.62 0.65 0.69 0.72 0.75 0.78 



crossover scale is given by T* = ^C{t_L/t)^^^^~^\ with 
e = {Kp + l/Kp)/A - 1/2 (~ 0.03 here). The interac- 
tions reduce its value compared to the naive one tj^/n. 
We cannot meaningfully test this renormalization effect 
because of the small values of 9 in the Hubbard model. 
Our results are consistent with T* ~ Ct^/n with a pro- 
portionality constant C ~ 0.5. In the low temperature 
FL regime, we find that the location of the FS is es- 
sentially unaffected by interactions. The quasi-particle 
residue Z is not uniform along the FS. Its dependence 
on fc^ is displayed in Table fl Z is larger for k±^ close to 
TT, while it is somewhat smaller for small fc_L. This is not 
a very big effect however, and it would be inappropriate 
to speak of "hot spots" Q in that region (even though 
it may be indicative that such a phenomenon appears 
for stronger interactions). Note that the single-chain ap- 
proximation E = TiiD used by many authors would lead 
to a much stronger variation of Z along the FS. Our re- 
sults, in agreement with those of ||l2[, show that Z is 
more uniform than this approximation would suggest. 

We now turn to the commensurate (half-filled) case, 
and present results for U = 0.65 M^. By fitting the self- 
energy to an analytical form, we first checked that our 
numerical results are consistent with a Mott insulating 
state with a gap /^id/W = 0.1, in agreement with exact 
results. Turning on t^, we estimate the i^-dependence of 
the gap by performing simulations at one of the lowest ac- 
cessible temperatures T/W = 1/40. The gap vanishes for 
t±/W ~ 0.07. Hence our method captures the insulator- 
to- metal transition induced by transverse hopping. The 
effective Kp, shown in Table pi is also a clear indicator 
of the MI {Kp = 0) to FL {Kp — 1) transition, in agree- 



TABLE IL Effective Kp at half-filling, as a function of t±/W 
for U/W = 0.65 and T/W = 1/40. 

t±/W 0.00 0.04 0.07 0.11 0.14 0.16 0.18 
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FIG. 3: Inter-chain optical conductivity at half-filling for 
U = 0.65W, P = AQ/W and ix = 0.14 and 0.07 (solid and 
dashed lines) and for U = l.QW , (3 = 40/1^ and tx = 0. 
(dotted line). 



ment with the qualitative expectations of Fig. nl At this 
low temperature, the intermediate LL regime is too nar- 
row to be visible. The location of the transition is in 
reasonable agreement with the naive criterion Aic ^ t*^ 
P, M, O] . Fig. ^ displays the inter-chain optical conduc- 
tivity for several values of t±. In the MI phase, a±{uj) 
shows a gap, followed by an onset of absorption starting 
at approximately the gap and extending up to a scale 
of order U , where a broad second peak is apparent. A 
low-frequency Drude peak develops as the insulator to 
metal transition is crossed. Close to the transition, the 
weight in the Drude peak is small, while the Hubbard 
band feature is still clearly visible and carries a signifi- 
cant part of the spectral weight. The results for the NMR 
relaxation rate 1/Ti are displayed in Fig. U for a value 
of t^/W — 0.11. At low temperatures the Korringa law 
\/{TiT) — const is recovered, while at higher temper- 
atures 1/TiT is temperature dependent. The observed 
T-dependence is consistent with l/{TiT) ~ T^"-^ |] 
with a T-dependent exponent Kp < 1, indicative of a 
LL evolving gradually into a FL (cf. Fig. ||). As for the 
doped case, the FS is nearly indistinguishable from the 
non-interacting case (see Fig. H). The QP residue Zk^^ 
depends only very weakly on the location on the Fermi 
surface as shown in Table pH. The data may suggest very 
shallow minima, which in contrast to the doped case, are 
found at k± '^ ±tt/2, corresponding to a vanishing inter- 
chain kinetic energy. We would like to contrast our find- 
ings with the conclusions drawn from the "single-chain" 
(RPA) approximation |p|, p3: S = Sid. As recently 
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FIG. 4: l/(riT) versus T in the half-filled case with ix/W = 
0.11, U/W = 0.65. The crossover from a FL at low-T to a LL 
for T /W > 0.025 is clearly apparent. 



TABLE III: QP weights Z{k^) for different points on the FS 
(half-filled case, tj_ = 0.14W, U/W = 0.65). 

fcx/TT 0.23 0.38 0.50 0.62 0.77 
Z{ki_) 0.79 0.77 0.76 0.77 0.79 



pointed out |2C(], this approximation predicts that the 
FS close to the metal-insulator transition should consist 
of disconnected electron and hole "pockets" , as depicted 
schematically on Fig. ^. This is because the self-energy 
of a MI diverges at low-frequency for k = 7r/2, so that 
no FS point can correspond to this value of A;. In con- 
trast, our data find a conventional connected FS down 
to the transition point (up to our numerical accuracy). 
This shows that it is crucial, as done in ch-DMFT , to 
take into account the feedback effects of the interchain 
hopping in the self-energy. 

In conclusion, we have shown in this paper that the 
ch-DMFT approach is a tool of choice for the study of 
the crossovers and insulator to metal transitions induced 
by transverse hopping in quasi onc-dimcnsional systems. 
Several features found in our study are qualitatively rem- 
iniscent of experimental observations on organic conduc- 
tors, such as: the deconfinement transition itself, the 
crossover from a LL at high-T to a FL at low-T (as 
revealed e.g in NMR) and most significantly the coex- 
istence of a Drude feature with small spectral weight and 
of a Hubbard-band feature in optical conductivity. Ex- 
tensions of the model including on-site and nearest neigh- 
bor interactions would allow to study the deconfinement 
transition at quarter-filling and thus to make a more real- 
istic ||3|, |9) comparison with experimental results. Though 
this is computationally much more demanding, we hope 
to use ch-DMFT in this context in future work. 
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FIG. 5: FS in the half-filled case with t^/W = 0.14, U/W = 
0.65 (circles), compared to the FS of the non-interacting case 
(dotted line) and of the purely Id case {t± — -dashed-.). The 
solid line depicts schematically the Fermi surface obtained 
when making the approximation E — T,io |20|. 
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